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ABSTRACT 

We present a pre-discovery iJ-band image of the HR 8799 planetary sys- 
tem that reveals all three planets in August 2007. The data were obtained 
with the Keck adaptive optics system, using angular differential imaging and 
a coronagraph. We confirm the physical association of all three planets, includ- 
ing HR 8799d, which had only been detected in 2008 images taken two months 
apart, and whose association with HR 8799 was least secure until now. We con- 
firm that the planets are 2-3 mag fainter than field brown dwarfs of comparable 
near-infrared colors. We note that similar under-luminosity is characteristic of 
young substellar objects at the L/T spectral type transition, and is likely due to 
enhanced dust content and non-equilibrium CO/CH4 chemistry in their atmo- 
spheres. Finally, we place an upper limit of >18 mag per square arc second on 
the >120 AU if-band dust-scattered light from the HR 8799 debris disk. The 
upper limit on the integrated scattered light flux is 10~ 4 times the photospheric 
level, 24 times fainter than the debris ring around HR 4796A. 



Subject headings: stars: low-mass, brown dwarfs — planetary systems — stars: in- 
dividual (HR 8799) — instrumentation: adaptive optics 
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INTRODUCTION 



The three-planet system around the A5V star HR 8799 (IMarois et al.ll2008bl . henceforth, 
M08) is the first directly imaged extrasol ar multi-planet sy stem. Along with a reported 
extrasolar planet around Fomalhaut (A4V; iKalas et al.l 120081 ) . HR 8799b, c, and d are also 
the first bona-fide extrasolar planets directly imaged around stars other than the Sun. The 
HR 8799 and Fomalhaut planetary systems share two important similarities. First, in both 
cases the planets are widely separated from their host stars, at projected separation s between 
24-1 1 9 AU. And second, both h osts are A stars surrounded by cold debris disks ( lAumann 
19851 ; ISadakane fc Nishidaj 1 1 98 6t ) . circumscribing the planetary systems. 



one of the most massive 


detected by IRAS and 


2.3 x 10" 4 ; 


Moor et al. 




2006; 


Rhee et al. 


2007) 



8 x 10 5 ; iRhee et al.l 120071 and references therein). 



is a factor of several brighter (Lir/L* 
than that around Fomalhaut (Lm/I, = 
While Fomalhaut 's disk has been spati ally resolved at wavelengths ranging fr om the visible 
fcalas et al.ll2005h to the submillimeter ^Holland et aljll998jjMarsh et alEooi l. the HR 8799 
debris disk remained unresolved until recently ( Su et al. 20091 . henceforth, S09). 



Prior to the discovery of the planet s around HR 8799 by M 08, we targeted the star with 
the Keck adaptive optics (AO) system (jWizinowich et al.l 120001 ) to detect the debris disk in 
scattered li ght. The strength of its IR e xcess and a detection in the submillimeter («10 mJy 
at 850 /im; I Williams fc Andrews! 120061 ) indicated a substantial optical depth, which given 
suitable disk viewing geometry could result in a detection of dust-scattered light. The data 
presented here did not produce the sought-after scattered light detection, but have revealed 
all three known planets around HR 8799. In particular, we report the earliest-epoch image 
of the closest-in planet, HR 8799d. Prior epoch de tectio ns of HR 8799b a nd/or HR 8799c 
have been published in M08. iFukagawa et al.l (120091 ) . and lLafreniere et al.l (120091 ). 



OBSERVATIONS 



We observed HR 8799 (H = 5.28 mag; 39.4 pc) with the Keck II AO system on 2 August 
2007. Atmospheric seeing was variable, ranging from 0'/6 — 1'.'4 in the visible. Cloud cover was 
pervasive and ranged from contributing ~1 mag in visual extinction to at times completely 
preventing tele scope operations. O bservations were taken at H band in angular differential 
imaging (ADI; IMarois et al.l 120061 ) mode, by fixing the telescope pupil rotator to optimize 
contrast and the detectability of non-face-on disks. We used the NIRC2 near-IR camera in 
its 10 mas pix -1 ("narrow") scale for optimal spatial sampling with the ADI technique. We 
targeted HR 8799 for a total of five hours around transit with frequent weather interruptions. 
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We obtained 114 one-minute exposures, of which 73 (i.e., 1.22 hrs total) were adequate 
for high-contrast science. All were taken with the star covered by the l'.'0-diameter semi- 
transparent NIRC2 coronagraph, and none were saturated. The LARGEHEX hexagonal 
pupil plane stop within NIRC2 was used to mask the edges of the primary mirror. Several 
unocculted AO-corrected images (5 ms, 1000 coadds per integration) were acquired during 
periods of stable seeing to calibrate the coronagraph's transmission. At flux levels <8,500 
cts pix -1 coadd -1 HR 8799 was in the linear response regime of the camera in all unocculted 
exposures. 

Sets of five 60 s coronagraphic exposures of a point-spread function (PSF) reference 
star, HR 8788 (F6V, H = 5.13 mag), were obtained at three different times during the 
night for additional PSF calibration. Images of the dome-illuminated interior with the l'.'O 
coronagraph in place were used for flat-field calibration. One-minute dark frames were taken 
at the beginning of the night. 



3. DATA REDUCTION 

The ADI sequences of exposures were dark-subtracted, flat-fielded, bad-pixel corrected, 
and reduced using cus tom-made ADI and LOCI ("locally optimized combination of images"; 



Lafreniere et al.ll2007l ) routines. The images were then spatially cross-registered and averaged 
to produce the final image shown in Figure [TJ Precise frame registration was possible because 
of the partially transmissive property of the NIRC2 coronagraph, which allowed alignment 
on HR 8799 itself. The final image is comparable in contrast to those published in M08, and 
all three known planets are detected. 

ADI observations optimize contrast for non-azimuthally symmetric circumstellar struc- 
ture, such as point source companions or a non-face-on disk. However, because the HR 8799bcd 
planetary system is seen near pole-on, and given the likelihood of planar alignment between 
the planetary orbits and the debris disk, we also independently reduced our data using 
the PSF reference star exposures. To mitigate the effect of variable AO performance we 
applied a radially-dep endent scaling to the reference PSF, enhancing a procedure used in 



Metchev et al.l (120051 ). More precisely, after registering the science and the reference-star 
exposures, we fit a linear relation fij(r) to the ratio %(r) = Si(r)/pj(r) of the radial profiles 
of each science exposure i with the nearest PSF exposure j. We then scaled the PSF by 
fij(r), and subtracted the thus re-normalized PSF from the science exposures before aver- 
aging. The linear fit to %(r) was obtained at 5r = 80 mas steps over overlapping 3"-wide 
search annuli. The PSF subtraction was thus aimed at revealing high order (quadratic or 
more) radial structure with widths up to ~3". The PSF-subtracted exposures were rotated, 
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Fig. 1. — Keck AO iJ-band image of the HR 8799bcd planetary system from 2 August 2007. 
The identification of HR 8799d in this image is the earliest by almost a year. 



- 5 - 



aligned, and then averaged and 3cr-clipped to produce the final image. 

Even with the use of a separate PSF reference star, we did not detect scattered light 
from the debris disk around HR 8799. This type of data reduction does show HR 8799b and 
c, although it does not reveal the innermost planet d. 



4. ANALYSIS AND DISCUSSION 

4.1. Astrometry of the HR 8799 Planetary System 

We measured the positions of the three planets by fitting two-dimenational Gaussians to 
their profiles in the reduced ADI image. Their relative astrometry with respect to HR 8799 is 
listed in Table [TJ We adopted a pixel scale of 9.963 ±0.005 mas pix -1 for the NIRC2 narrow 
camera, a 0.13 ± 0.02° o ffset from North for the camera columns, and corrected for the 



known camera distortion (IGhez et al.l 120081 ) . Astrometric errors were derived by calculating 



the registration bias from the residual image noise found near the location of each planet. 

The measured positions confirm the associations of all three planets with HR 8799, and 
are consistent with counter-clockwise orbital motion along nearly face-on orbits. Due to a 
merely two-month time span between the discovery and confirmation epochs, HR 8799d had 
the lowest significance (6<x) of physical association with the host star in M08. It is now con- 
firmed as a bound object at the ~13<7 level from the ensemble of astrometric measurements 
from four independent epochs (Fig. [2]). HR8799d was also the planet with the least accurate 
orbital motion estimate in M08 (42 ± 27 mas yr _1 ). The longer time-line afforded by the 
present observations indicates an orbital motion of 41 ± 9 mas yr _1 (1.6 ± 0.4 AU yr _1 ). For 
a semi-major axis of 24.3 AU and a primary mass of 1.5 M , the orbital period is 98 years 
and the expected motion along a circular face-on orbit is 1.57 AU yr -1 . 



4.2. Photometry of HR 8799bcd 



Photometry on the planets was performed on the reduced ADI image as described in 
M08. We calibrated the planet fluxes relative to that of HR 8799 (H = 5.280 ± 0.018 mag; 



Skrutskie et al.ll2006l ) seen through the l'.'O-diameter NIRC2 coronagraph, for which we mea- 
sured an attenuation of AH = 7.78 ± 0.10 mag from compariso ns with the short unoc- 
culted exposures: in agreement with earlier calibrations performed in lMetchev &: Hillenbrand 
( 120091 ). The ADI/LOCI algorithm and field rotation during exposures were also calibrated by 
introducing artificial field-rotated planets using the unocculted PSFs. The apparent i/-band 
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Fig. 2.— Orbital motion of HR 8799d between 2 August 2007 (position 1) and 12 September 
2008 (position 4) relative to the position of HR 8799. The dotted line with perpendicular 
dashes at the corresponding observational epochs traces the expected change in position of 
HR 8799d due to the proper and parallactic motion of the primary if the planet were an 
unrelated background object. 
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magnitudes for HR 8799bcd are shown alongside those from M08 in Table [TJ 



Our photometry for HR 8799b and c is in good agreement with that published in M08. 
The slightly better precision of our measurements is due to the simultaneous calibration with 
respect to HR 8799 itself. M08 did not use a coronagraph for their H-band imaging and 
calibrated their photometry from non- simultaneous data: from short sequences of shallow 
unsaturated exposures taken between the long sequences of saturated deep exposures. Our 
setup avoids non-simultaneous calibration of the data, which can result in difficult to quantify 
seeing/sky throughput variations. HR 8799d is 0.4 mag (1.4<r) fainter than reported in M08. 
While the difference could be random or could indicate real variability in HR 8799d, it could 
also be caused by a possible decrease in throughput of the optical system in close proximity 
(~0'.'l = 10 pix ftj 3 FWHM) to the coronagraph. if-band photometry of HR 8799b obtained 
at earlier epochs (see Table [T]) is consistent with the measurements presented both here and 
in M08. Hence, we do not see evidence for significant variability in any of the planets. 

Compared to (mid-T) brown dwarfs of similar intrinsic brightness (Mh ~ 15 mag), 
HR 8799b, c, and d appear unusually redder, by > 1 mag in H — K, and are matched onl y 
by the ~8 Myr-old 2MASS 1207334-393254B (H - K ~ 1.2 mag; Ichauvin et aD l2004h . 
Conversely, all four of these planetary-mass objects are 2-3 mag fainter at H band than 
comparably red >1 Gyr-old mid- to late-L field dwarfs. 

Similar under-luminosity with respect to the field substellar population, although not 
as extreme (by ~ 1 mag at H), is already known in ~300 Myr-old brown dwarfs at the 
L/T spectral type transition, wh ere it has been interpreted as evidence of unusual ly low 
effective temperatures (<1200 K; iMetchev fc Hillenbrand! 120061 ; iLuhman et al.l 120071 ). The 
warmer than expected spectroscopic appearance of such young cool substella r objects may be 



linke d to unusually high dust content in their low-gravity atmospheres (M08; lLafreniere et al. 



20091 ). In addition, the fact that the phenomenon is observed in all young L/T-transition 
dwarfs and only in them (and, to a lesser extent, in a moderately young mid-T brown dwarf 



2008^ hints at an augmented role of CO/CH4 non-equillibrium chemistry 
( ISaumon et al.l 120061 ) at low substellar surface gravities. This hypothe sis is further sup - 



binary; iLiu et al. 



ported by the recent non-detection of the HR 8799 planets at M-band by iHinz et al.l (120091 ) 
suggesting unusual strength of the CO fundamental bandhead at 4.7/xm. 



4.3. Upper Limits on Outer Planets 



The 10"2 x 10'.'2 NIRC2 field of view (FOV) and a ~1" off-center positioning of HR 8799 
allows us to search for planets out to ~4" from the star. We do so by considering the pixel 
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Fig. 3. — Detection limits for point sources (in magnitude contrast AH with respect to the 
primary; solid line) and extended emission (in magnitudes \in per square arc second; dashed 
line) around HR 8799 (H = 5.28 mag). The three detected planets are marked. The point 
source limits are from the ADI image reduction approach. The extended source limits are 
obtained with the use of HR 8788 for PSF calibration. 
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flux scat ter in our AP I -reduc ed final image as a function of radial separation from the star 
(Fig. [3]). iMarois et al.l (l2008al ) have shown that the residuals obtained from averaging many 
PSF-subtracted exposures taken with the ADI technique behave in a Gaussian fashion, and 
so the detectability of faint point sources in the contrast-limited regime can be rigorously 
determined from the scatter of the pixel fluxes. To account for the spatial correlation among 
neighboring pixels on the scale of a PSF, we convolve the ADI-reduced image with the PSF 
(FWHM = 3.4 pix = 34 mas at H band after using a high-pass filter in the data reduction). 

We adopt a 4<r detection limit for point sources as a l evel of significance that closel y 
matches the ability to visually identify faint companions (IMetchev &: Hillenbrand! 120091 ) . 
Incidentally, we note that 4<r is the approximate level at which we detect HR 8799d (Fig. [3]). 
We do not find any H < 22 mag planets beyond th e 68 AU projected separation of HR 8799b, 
complementing results from I Close fc Males! (120091 ) . who find no H < 22 mag planets between 
5-15" (200-600 AU). There fore, for an adopted 60 Myr age of HR 8799, we can exclude 
additional >3 M Jup planets (IBurrows et al.ll2001uBaraffe et al.ll2003l ) out to at least 600 AU. 



4.4. Debris Disk 

Given the expectation for a face-on debris disk around HR 8799, we only consider the 
PSF subtraction performed with the use of a separate PSF star, as that is better suited for 
revealing uniform extended emission. We disregard the inner region of the planetary system, 
and focus only on angular separations >1"5 (>60 AU), where the ratios qij(r) of the HR 8799 
to PSF radial profiles behave in a nearly linear fashion out to the edge of the FOV. S09 have 
recently resolved extended 24 /^m and 70 fim emission from dust orbiting HR 8799 at such 
separations, which they model as arising from a cool planetesimal belt beyond 90 AU. 

Analogous to the procedure adopted for point source detection (§ 14.31) . we convolve our 
reduced image with the PSF and impose a 4cr detection threshold. We add in quadrature 
to this (random) source of error an estimate of the systematic uncertainty cr sys (r) in fitting 
the large-scale wings of the PSF under varying seeing conditions, i.e., in determining ftj(r). 
a sys (r) is obtained as the standard deviation of the mean of fij(r)pj(r), and we find that it 
contributes approximately 50% of the error budget. 

We do not detect any if-band scattered light from the HR 8799 debris disk down to 
a limit of /i# «18 mag (65 //Jy) per square arc second at 3-6" (120-240 AU; Fig. [3]). In 
particular, we do not detect near-IR scattered light from the resolved far-IR disk seen in S09, 
although it is possible that our 3"-wide search annulus may be over-resolving the observed 
extended emission. Given our 10"2 x 10"2 FOV, a wider search annulus was not practical. 
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It is instructive to compare the upper limit on the surface brightness of the HR 8799 
debris disk to other debris disks seen in scattered light, and in particula r to th e ring around 
HR 4796A, resolved at 1.6 fim (F160W) with the HST bv lSchneider et al.Ul999h . For simplic- 
ity of comparison, we assume that the HR 8799 debris disk also has a ring-like morphology, 
with an aspect ra tio (width : radius ~ 1:4) similar to that of HR 4796A and Fomalhaut 
(IKalas et al.ll2005l ). and to a model for the cool planetesimal belt shown in Figure 9a of S09. 
Given the inferred near-pole-on orientation of the planetary orbits (M08) and the circular 
symmetry of the resolved far-IR emission (S09), we assume that the HR 8799 debris disk is 
viewed nearly pole-on. 

For a hypothetical disk inner edge at either 89 AU, in ~3:2 resonance with HR 8799b 
and coincident with the 90 AU inner edge for the cool planetesimal belt from S09, the 4<r 
upper limit on the integrated if -band disk surface brightness is 15.43 mag (0.69 mJy). That 
is, the fraction of dust-scattered li ght is / sra ,t < 1 0~ 4 : a factor of >24 fainter than for the 



HR 4796A disk (f BCat = 2.4 x 10~ 3 ; ISchneider et al.lll999h . If the HR 87 99 debris disk were 



more extended, similar to that around HD 107146 (width : radius ~ l:2; lArdila et al.ll2004l ) 
or as in the preferred model for the cool planetesimal belt in S09 (see their Figure 9b), then 
/ scat < 2 x 10- 4 . 



Given the strength of the HR 8799 IR excess (f m = 2.3 x 10" 4 ), the ratio of the disk's 
IR to scattered lig ht is fm./ /scat > 2.3. This compares to fm/ fsc&t = 1-8 for HR 4796A 



(/ m = 4.4 x 10~ 3 ; IRhee et al.ll2007h . The higher f m /f Bcat ratio for HR 8799 marginally 



excludes the presence of efficient near-IR scatterers between 90 and 240 AU, consistent with: 
(1) the minimum grain size of several microns set by radiation pressure (in the assumed 
absence of gas), (2) S09's deduced absence of <10 fim grains in the cool pl anetesimal belt, 



and (3) the inferred lack of sub-micron grains in the disk around HR 4796A (ISchneider et al. 
1999h . 



5. CONCLUSIONS 



We have presented an year-2007 image of all three known HR 8799 plan ets, including 



a detec tion of the inner-most planet HR 8799d that precedes its discovery in iMarois et al. 



(j2008bl ) by one year. Our data exclude the presence of additional >3 M Jup outer planets 
between 68 and 160 AU from the star. We do not detect scattered light from the debris 
disk, and place an upper limit of ~18 mag arcsec -2 on its if -band surface brightness at 3-6" 
(120-240 AU) f rom the star. Based o n a more robust if-band photometric calibration than 
was possible in IMarois et ajj (l2008bl ). we confirm that HR 8799b, c, and d are very under- 
luminous or redder compared to field brown dwarfs of similar colors or intrinsic luminosities. 
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We note that this is a distinctive feature of all young substellar objects near the L/T transi- 
tion, and speculate that it may be due to enhanced dust content and a significant departure 
from CO/CH4 chemical equilibrium in cool low-surface gravity substellar atmospheres. 
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Table 1. Astrometry and iJ-band Photometry for HR 8799bcd from 2 August 2007 





Aa 


A<5 


H 




Planet 


(arcsec) 


(arcsec) 


(mag) 


(mag) 


b 


1.522 ±0.003 


0.815 ±0.003 


18.06 ± 0.13 


17.85 ±0.17 


c 


-0.672 ± 0.005 


0.674 ± 0.005 


16.95 ± 0.16 


16.91 ±0.17 


d 


-0.170 ±0.008 


-0.589 ± 0.008 


17.25 ± 0.20* 


16.84 ±0.22 



t Apparent H-band magnitudes from M08. 

tThe magnitude of HR 8799d may carry additional systematic uncertainty 
because of its proximity to the coronagraphic mask. 



Note. — Additional prior-epoch _ff-band phot ometry for HR 8799b: F160W 
= 18.54 ± 0.12 mag (fro m the HST in 1998; [ Lafreniere et a.1.1 [iooih ; H = 
18.07 ± 0.25 mag (in 2002; iFukagawa et al-lbooal 



